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Open 



Although all plant and animal species harbor microbial symbionts, we know surprisingly little about 
the specificity of microbial communities to their hosts. Few studies have compared the microbiomes 
of different species of animals, and fewer still have examined animals in the wild. We sampled four 
pond habitats in Colorado, USA, where multiple amphibian species were present. In total, 
32 amphibian individuals were sampled from three different species including northern leopard 
frogs (Lithobates pipiens), western chorus frogs (Pseudacris triseriata) and tiger salamanders 
(Ambystoma tigrinum). We compared the diversity and composition of the bacterial communities on 
the skin of the collected individuals via barcoded pyrosequencing of the 16S rRNA gene. Dominant 
bacterial phyla included Acidobacteria, Actinobacteria, Bacteriodetes, Cyanobacteria, Firmicutes 
and Proteobacteria. In total, we found members of 18 bacterial phyla, comparable to the taxonomic 
diversity typically found on human skin. Levels of bacterial diversity varied strongly across species: 
L. pipiens had the highest diversity; A. tigrinum the lowest. Host species was a highly significant 
predictor of bacterial community similarity, and co-habitation within the same pond was not 
significant, highlighting that the skin-associated bacterial communities do not simply reflect those 
bacterial communities found in their surrounding environments. Innate species differences thus 
appear to regulate the structure of skin bacterial communities on amphibians. In light of recent 
discoveries that some bacteria on amphibian skin have antifungal activity, our finding suggests that 
host-specific bacteria may have a role in the species-specific resistance to fungal pathogens. 
The ISME Journal (2012) 6, 588-596; doi:10.1038/ismej.201 1 .129; published online 29 September 2011 
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Introduction 

All species of plants and animals harbor assem- 
blages of microbes and these symbionts often have 
vital roles in the health of their hosts. However, 
there have been few comprehensive studies asses- 
sing whether different species harbor distinct 
microbial symbionts. Work on mammals (Ley 
et al, 2008), primates (Yildirim et al, 2010) and 
plants (Redford et al., 2010) suggests that there is 
some degree of host specificity of microbial sym- 
bionts across taxa. However, such studies have not 
yet been performed in wild populations, nor have 
these studies been designed to distinguish environ- 
mental factors from innate host factors. For different 
species in the wild, the question of host-specific 
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microbial assemblages becomes difficult to address 
because species often occupy different habitat 
niches, and different environmental conditions 
may confound questions of host specificity. For 
several reasons, we have identified pond-dwelling 
larval amphibian communities as a model system to 
address whether different, co-habiting species har- 
bor unique microbial communities. First, popula- 
tions of larval amphibians observed before 
metamorphosis (for example, tadpoles), develop 
within a single pond and therefore have a known 
environmental origin. Second, ponds can harbor 
several species of co-habiting larval amphibians that 
occupy the same environmental niche during their 
development. Third, different ponds that harbor 
multiple species of larval amphibians can serve as 
replicates to examine patterns of host-associated 
microbial community assembly across different 
sites. 

We focused specifically on skin-associated bacter- 
ial communities of amphibians as these commu- 
nities likely provide the first line of defense against 
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pathogens and, for co-habiting species, are exposed 
to the same environmental sources of microbes. One 
specific skin pathogen of global significance is the 
fungus Batrachochytrium dendrobatidis [Bd), which 
has a key role in amphibian declines worldwide 
(Skerratt et al., 2007; Fisher et al, 2009). Amphibian 
species, and some populations within species, vary 
in their susceptibility to disease caused by this 
fungal pathogen (Kilpatrick et al., 2010). While there 
are many factors that may have a role in amphibian 
disease resistance including host genetics, immu- 
nology (Carey et al., 1999), skin peptides (Rollins- 
Smith 2009) and environmental factors (Lips et al., 
2008), microbial communities living on the skin of 
amphibians may also have a crucial role in host 
defense mechanisms. Certain bacteria isolated from 
amphibian skin are known to produce fungicides 
and can inhibit growth of Bd (Harris et al., 2006) and 
even reduce amphibian mortality under controlled 
experimental conditions (Harris et al., 2009b). 
However, it is not known if different species of wild 
amphibians can select unique skin microbial assem- 
blages from exposure to the same microbial back- 
ground. If specific microbial assemblages are 
involved in resistance to pathogens (such as Bd), 
and amphibian species vary in their susceptibility to 
disease, then we would expect that unique skin 
microbial composition across species or populations 
might provide a mechanism underlying this 
observed pattern of differences in susceptibility. 
Our goal in the present study was to use co-habiting, 
pond-dwelling larval amphibians to determine 
whether different species of amphibians harbor 
unique communities of bacteria on their skin. If 
they do, it suggests that species-specific bacterial 
communities have the potential to contribute to 
differences in the resistance of individual species to 
skin pathogens. 

Teasing apart how symbiotic microbial commu- 
nities are influenced by factors such as the environ- 
ment, diet or innate host traits can now be addressed 
with the use of newer high-throughput survey 
techniques such as barcoded pyrosequencing of 
the 16S rRNA gene. These methods allow for 
observing the composition and relative abundance 
of microbes with greater accuracy than previous 
methods used to estimate bacterial community 
structure (Costello et al., 2009). By applying these 
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newer techniques in pond-dwelling amphibians, we 
hope to gain insight on the host-specific nature 
of skin microbes independent of environmental 
confounds as well as build a foundation to under- 
stand the potential for these microbes to have a role 
in host resistance against pathogens. 



Materials and methods 

Amphibian species and field collections 
Four pond habitats in Colorado were identified in 
2008 that each harbored at least two amphibian 
species that are common in Colorado and co-habit in 
ponds during the summer breeding season. Amphi- 
bian species sampled included Lithobates [ = Rana) 
pipiens (northern leopard frog), Pseudacris triseriata 
(western chorus frog) and Ambystoma tigrinum 
(tiger salamander). Each pond was less than one 
hectare in surface area and located in a separate 
watershed. Site characteristics and the number of 
individual amphibians sampled per species per 
pond are described in Table 1. Amphibians were 
sampled for skin bacteria during August 2008 and 
all individuals from a given pond were sampled on 
the same day. Pre-metamorphic, larval amphibians 
were targeted because they develop over one season 
within a single pond, therefore, we can ensure that 
our cross-species comparisons only include indivi- 
duals that shared the same environment and micro- 
bial inocula during their development. All three 
species named above are pond-breeders that lay eggs 
in early summer, which hatch larval swimming 
stages that develop over the course of the summer 
into metamorphic stages (Hammerson, 1999). This 
system allows us to control for confounding envir- 
onmental conditions and isolate the effect of host 
species. 

Larval amphibians were captured from the water 
using hand-held dipnets or a blocking net seine. 
Within a given species, larval individuals of a 
similar size were selected for skin microbe sam- 
pling. Larval amphibians were held temporarily 
(<10min) in plastic containers with pond water 
before sampling. In total, 32 individual amphibians 
were sampled (7L. pipiens, 14 P. triseriata and 12 
A. tigrinum; see Table 1) from four ponds. Fresh latex 
gloves were used for every individual amphibian 



Table 1 Amphibian sampling scheme 



Pond 


# Leopard frogs 


# Chorus frogs 


# Tiger salamanders 


Pond elevation 


Pond location in 


Latitude longitude 


name 


sampled 


sampled 


sampled 


(m) 


Colorado 


(UTM) 


Monster 


4 


2 


4 


2357 


Meeker 


40.010359-107.625487 


bull 














Chive 


0 


4 


4 


2598 


Nederland 


39.951877-105.463940 


Horse 


0 


4 


4 


1702 


Southwest Boulder 


39.946662-105.241441 


Eggleston3 


3 


4 


0 


1754 


Southeast Boulder 


39.927857-105.20637 



All three species co-occurred at one site and other sites had two species co-occurring. 
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handled. Before sampling, each individual was held 
and rinsed with 100 ml of sterile water three times to 
ensure that the skin sample primarily included skin- 
associated microbes rather than pond-associated 
material, including pond water, sediment and 
transient microbes (Culp et al., 2007; Lauer et al., 
2007). Prior studies by Culp et al. (2007) and Lauer 
et al. (2007) demonstrate that the composition of 
bacteria obtained from amphibian skin versus rinse 
water mostly differ and represent unique assem- 
blages, thus, we have some precedent to assume that 
most of the bacteria observed in this study are likely 
associated with the amphibians rather than transient 
bacteria from the environment. Immediately follow- 
ing rinsing, each amphibian was sampled using a 
sterile cotton-tipped swab brushed over the entire 
body of the amphibian for 30 s. The small size of 
larval amphibians (3-12 cm, depending on species) 
prevented us from swabbing a specific site on the 
body, so we opted to uniformly swab the entire body 
(head, mouth, tail, dorsal and ventral surfaces) in an 
effort to capture the total skin microbiome of each 
amphibian. All individuals were swabbed according 
to a standard protocol such that the swab was 
brushed in the direction from the head to the tail 
and the swab was turned slightly with each move to 
a new body surface (for example, dorsal stroke, 
rotate swab, then lateral body stroke) in order to 
capture each body surface on the swab without 
saturating the whole swab with one body part; this 
technique was manageable given the small size of 
the animals. Each swab was placed in a sterile vial 
and stored on dry ice. Within a few hours, each 
sample was then transferred to a — 20°C freezer for 
storage until DNA extraction. All animals were 
handled and released according to an approved 
Institutional Animal Care and Use Committee 
protocol. 



DNA extraction, bacterial 16S rRNA PCR amplification 
and pyrosequencing 

Bacterial DNA was extracted from each swab using 
the PowerSoil DNA isolation kit (MoBio Labora- 
tories, Carlsbad, CA, USA). For each of the 32 
samples, the swab was loaded into a bead tube from 
the extraction kit and the remaining procedures 
were performed according to the manufacturer's 
protocol. Bacterial community composition was 
determined using barcoded pyrosequencing, follow- 
ing the procedure described in detail in Fierer et al. 
(2008), which includes both the PCR conditions and 
error-correcting barcoded primer sequences. The 
primer set used was the 27F/388R that targets the 
V2 region of the 16S small subunit ribosomal gene. 
No template and template from blank filters were 
included as negative controls throughout the entire 
process from DNA extraction to PCR amplification 
to check for contamination. PCR amplicons were 
pooled from each sample to equimolar ratio in a 
single tube, and sequenced at the University 



of South Carolina Environmental Genomics Core 
Facility on a 454 Life Sciences Sequence GS FLX 
instrument (Roche, Florence, SC, USA). We obtai- 
ned an average of 1220 sequences per sample (range 
780-1510), with read lengths averaging 258 bp in 
length. Previous work has shown that reads of this 
length over the V2 region are sufficient for accurate 
taxonomic identification to at least the family level 
and for community characterization (Liu et al., 
2007). 

A portion of the extracted DNA was used to screen 
samples for Bd, the chytrid fungal pathogen. Frozen 
DNA was shipped to Pisces Molecular LLC (Boulder, 
CO, USA) for analysis using £?d-specific primers 
(Annis et al., 2004), PCR amplification and gel 
electrophoresis to examine for the presence or 
absence of bands. 

Sequence analysis 

All sequence analyses were conducted using 
the Quantitative Insights Into Microbial Ecology 
(QIIME) pipeline (Caporaso et al., 2010). In brief, the 
QIIME pipeline removed low-quality sequences 
(that is, those sequences < 200 bp in length, had a 
quality score <25, contained ambiguous characters, 
having an unreadable barcode or did not contain the 
primer sequence). Remaining sequences were clus- 
tered using CD-HIT (Li and Godzik, 2006), using a 
minimum coverage of 99% and minimum sequence 
identity of 97%. A representative sequence was 
chosen for each phylotype by selecting the longest 
sequence with the highest number of hits to 
other sequences in the same phylotype. The repre- 
sentative sequences were aligned using PyNAST 
(Caporaso et al., 2010) and taxonomy was assigned 
using the RDP classifier (Wang et al., 2007). A 
phylogenetic tree of all aligned sequences was cons- 
tructed with FastTree (Price et al., 2009) and used in 
all downstream phylogenetic community compar- 
isons. Community similarity (that is, (3-diversity) 
analyses were performed using the taxon-based 
Bray-Curtis dissimilarity metric, and the phylo- 
genetic UniFrac metric: the Bray-Curtis metric 
examines taxon overlap, whereas the UniFrac metric 
is a phylogenetic approach based on the fraction of 
shared branch length between two communities 
within a phylogenetic tree (Lozupone et al., 2007). 
Phylogenetic trees were constructed using FastTree 
(Price et al., 2009). All samples were rarified to 
750 sequences per sample to remove sample hetero- 
geneity, which impacts oc- and (3-diversity metrics 
(Lozupone et al., 2011). 

Statistical analyses 

To address whether bacterial phylotype richness 
varied among species and/or among pond sites, we 
conducted an analysis of covariance. Phylotype 
richness per frog served as the response variable 
(i7 = 32) and frog species and pond site served as 
two categorical covariables. This analysis was 
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performed using the statistical software JMP (SAS 
version 8, Cary, NC, USA). The data were normally 
distributed as determined by a Shapiro-Wilk W 
goodness-of-fit test. 

To determine the differences in bacterial commu- 
nity composition across individuals ((3-diversity), 
we used two different distance metrics: a Bray- 
Curtis similarity index (a taxonomic metric) and 
the UniFrac algorithm, which provides a measure 
of phylogenetic distance between communities 
from individual samples (Lozupone et al., 2007). 
P-Diversity patterns were visualized using a MDS 
(non-metric multidimensional scaling) ordination 
approach (Clarke and Warwick, 2001). To test for 
significant differences in community composi- 
tion, we used a two-way nested ANOSIM (analysis 
of similarity) to examine whether host amphibian 
species and/or pond site were significant predictors 
of variability in bacterial communities across 
individual amphibian samples (PRIMER-E). We ran 
two separate analyses, one with Bray-Curtis simi- 
larity as the response variable and one with 
unweighted UniFrac distances as the response 
variable. 

Lastly, to visualize the relative abundances of the 
most common bacterial groups across amphibian 
species per pond, we constructed a heat map figure. 
By removing all sequences that comprised ^5% of 
the total number of classified sequences across all 
samples, 19 phylotypes remained. By reducing the 
data set to the 19 most abundant phylotypes, 
visualization of the compositional differences was 
tractable and issues related to sequencing error are 
avoided (Quince et al., 2009; Huse et al., 2010; 
Kunin et al., 2010). 



Results 

From the rarefied 750 sequences per sample [n = 32), 
we found between 15 and 345 unique phylotypes 
per sample representing 18 different bacterial phyla. 
Diversity results thus apply only to the dominant 
members of the community, and it is possible that 
sequences scored as absent are in fact present at 
lower abundance; however, the communities are 
substantially different from one another when 
sampled at this depth. The dominant phyla across 
all samples were Acidobacteria, Actinobacteria, 
Bacteriodetes, Cyanobacteria, Firmicutes and Pro- 
teobacteria. A list of all phylotypes and their 
taxonomic affiliation is provided in Supplementary 
Information (Supplementary Table Si). 

Phylotype richness varied between the amphibian 
species (Figure 1). Northern leopard frogs had the 
highest bacterial phylotype richness relative to the 
first 750 sequences sampled (mean = 287 phylotypes 
per sample), western chorus frogs had a moderate 
richness (mean = 185) and tiger salamanders had 
the lowest (mean = 66). These phylotype richness 
patterns varied significantly among species 
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Figure 1 Phylotype richness per amphibian species. The 
number of unique phylotypes per individual amphibian, 
estimated from a sampling depth of 750 sequences per sample. 
Each symbol represents the estimate of unique phylotypes or 
operational taxonomic units (OTUs) per amphibian individuals at 
a given site. Different symbols correspond to different pond 
locations. 



(P<0.0001, r 2 = 0.60) but were not significantly 
different among ponds [P=0A7), as determined by 
an analysis of covariance. We cannot effectively 
estimate absolute levels of diversity because sequen- 
cing errors may lead to an inflation of phylotype 
numbers (Quince et al., 2009) and because deeper 
sequencing depths would undoubtedly increase the 
number of phylotypes identified. However, since all 
samples were compared at an equivalent sequencing 
depth, these results do suggest that there are strong 
differences in relative levels of bacterial diversity 
across the amphibian species sampled here. 

There were strong differences in the composition 
of the skin-associated bacterial communities found 
on the different amphibian species, a pattern 
evident in the ordination plot (Figure 2). These 
patterns are confirmed by the ANOSIM analyses; 
with both the taxon-based Bray-Curtis distance and 
the unweighted UniFrac phylogenetic distance 
metrics, host species identity was a highly signifi- 
cant predictor of skin bacterial community unique- 
ness [P< 0.001; Table 2). Host species explained a 
large proportion of the variance in the skin bacterial 
communities (86% using Bray-Curtis and 91% 
using unweighted UniFrac). Different pond sites 
were not significant predictors of skin bacterial 
communities on amphibians (Table 2), although 
variation among sites within species does exist, as 
evident from Figure 2. 

The dominant phylotypes across all amphi- 
bian samples and their relative abundances are 
presented in Figure 3, with further details provided 
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Figure 2 Ordination plot. |3-Diversity patterns were visualized using a MDS ordination approach with skin-associated bacterial 
community differences represented as Bray-Curtis distances. The 2D stress (or measure of distortion) in the MDS configuration was 
relatively low (0.13), so the 2D distance between points in the ordination plot is a good representation of the degree of similarity between 
each sample's bacterial community. Each point represents the skin bacterial community of an individual amphibian; color indicates 
species and shape indicates pond location. 



Table 2 Results of the two-way nested analysis of similarity (ANOSIM) 





Bray-Curtis similarity 


Unweighted UniFrac distance 




Global R 


Significance 


Global R 


Significance 


Between amphibian species 


0.868 


< 0.001 


0.914 


<0.001 


Between pond sites 


0.089 


0.391 


-0.178 


0.714 



We examined the differences among bacterial communities from different amphibian species nested within different pond sites using two 
different community similarity measurements. 



Amphibian species 



Tiger salamander 



Pond name 

Abundant Bacterial phylotypes: N (# frogs) 


Chive Pond 
4 


Horse Pond 
4 


Monster Bull 
4 


Acidobacteria-Acidobacteriales 




Actinobacteria-Actinomycetales 


0.03% 


0.07% 


1.33% 


Alphaproteobacteria-Caulobacteraceae 


0 


0 


0.10% 


Alphaproteobacteria-Sphingomonadaceae (Sphingomonas) 


2.21% 


0 


1.22% 


Bacteroidetes-Bacteroidaceae 


0.48% 


0 


2.15% 


Bacteroidetes-Flavobacteriaceae 


0.50% 


0.03% 


0.52% 


Bacteroidetes-Sphingobacteriales 


0.50% 


0.57% 


51.29% 


Betaproteobacteria-Comamonadaceae (Curvibacter) 


80.19% 97.14% 12.67% 


Betaproteobacteria-Hydrogenophilaceae 


0 


0 


0.17% 


Betaproteobacteria-lncertae sedis 5 {Ideonella) 


0.03% 


0 


0.10% 


Betaproteobacteria-Neisseriaceae 


° 


0 


0 


Cyanobacteria 


0.54% 


0.03% 


0.79% 


Deltaproteobacteria 


0.27% 


0.03% 


0.13% 


Firmicutes-Clostridiaceae 


2.79% 


0 


0.13% 


Firmicutes-Peptococcaceae 


2.86% 


0 


6.39% 


Gammaproteobacteria-Enterobacteriaceae 


0.27% 


0 


3.73% 


Gammaproteobacteria-Pasteurellaceae 


0 


0 


0 


Gammaproteobacteria-Pseudomonadaceae 


0.07% 


0 


0.10% 


Gammaproteobacteria-Xanthomonadaceae (Nevskia) 




0 


0.07% 


Total proportion of sequences represented in this subset: 


0.91 


0.98 


0.81 



Western Chorus frog 


Chive Pond 


Eggleston3 


Horse Pond 


Monster Bull 


4 


3 


4 


2 


2.48% 


0.46% 


3.67% 


0 


10.38% 


3.62% 


10.58% 


0.28% 


1.09% 


6.94% 


0.94% 


0 


7.41% 


8.86% 


3.97% 


4.11% 


0.33% 


0.04% 


0.96% 


0.14% 


3.10% 


1.23% 


0.97% 


1.43% 


1.83% 


0.50% 


0.96% 


0.14% 


12.02% 


42.38% 


20.46% 




4.11% 


2.17% 


0.62% 


0 


2.44% 


3.49% 


17.84% 


0 


0.11% 


0.54% 


0 


0 


4.69% 


1.79% 


1.81% 


0.14% 


0.59% 


0.32% 


0.35% 


0.35% 


0.11% 


0.05% 


0.15% 


0 


0.11% 


0 


0.08% 


0 


0.18% 


0.23% 


0.27% 


0 


0.52% 


4.68% 


0.03% 


0 


7.06% 


2.58% 


1.20% 


0.07% 


8.08% 


0.55% 


7.56% 


0 


0.67 


0.80 


0.72 


0.98 



Northern Leopard frog 



Eggleston3 
3 



Monster Bull 
4 



3.58% 


1.16% 


8.16% 


5.66% 


0.68% 


0.43% 


5.88% 


3.97% 


2.44% 


0.12% 


0.92% 


3.65% 


3.25% 


1.73% 


13.73% 


27.56% 


1.04% 


0.91% 


2.78% 


1.73% 


0.21% 


3.85% 


5.74% 


7.86% 


3.26% 


2.19% 


0.10% 


0.28% 


0.34% 


0.63% 


0.10% 


0 


0.24% 


0 


1.02% 


3.32% 


1.64% 


1.62% 



Figure 3 Heat map. A color-scale heat map demonstrates the relative abundance of the 19 most common bacterial phylotypes across 
species and ponds; red cells indicate higher proportional abundance and blue cells indicate lower proportional abundance. The relative 
abundance of each phylotype was converted to a proportion of the total number of classified sequences. The total number of unclassified 
sequences summed together for each amphibian sample comprised on average only 5.7% ( ± s.d. = 0.06) of the sequences per sample and 
were removed from the table. In order to condense the bacterial groups and present the most common phylotypes in the heat map, we 
further removed all phylotypes that comprised ^5% of the total number of classified sequences across all samples. The proportion of 
sequences that are represented in the heat map for each species in each pond site varied between 55-98%. Nineteen unique phylotypes 
remained and the average proportion was calculated per amphibian species per pond. 
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in Supplementary Table Si. The most dominant 
phylotype across all three amphibian species was 
classified as belonging to the genus Curvibacter 
[Betaproteobacteria: Comamonadaceae), which 
accounted for 12-97% of the sequences from 
individual samples, depending on the species and 
pond location (Figure 3). A Bacteroidetes phylotype 
belonging to the group Sphingobacteriales was 
dominant on tiger salamanders from one pond, 
comprising 51% of the community, but otherwise, 
the remaining phylotypes observed on salamanders 
comprised ^6% of the bacterial sequences. 
A species of Ideonella [Betaproteobacteria] was 
abundant on chorus frogs from one pond, compris- 
ing 18% of the community and another phylotype 
belonging to the Actinobacteria [Actinomycetales) 
was also abundant on chorus frogs at two of the 
pond sites. Leopard frogs had many phylotypes 
present at lower abundances. For instance, several 
phylotypes comprised between 3% and 8% of the 
leopard frog skin communities, and these phylo- 
types belonged to a range of different phyla and 
subphyla. 

One individual western chorus frog [P. triseriata) 
from Chive Pond was found to be positive for Bd, the 
chytrid fungal skin pathogen. Due to a lack of other 
positive individuals in this data set, we cannot draw 
any conclusions about the possible interactions of 
chytrid infection with the skin bacterial commu- 
nities. 

Discussion 

Our understanding of the microbes associated with 
plants and animals, including humans, is in a phase 
of rapid discovery owing, in part, to the develop- 
ment of high-throughput sequence-based appro- 
aches that allow us to document the structure and 
composition of numerous complex bacterial com- 
munities simultaneously. With improved commu- 
nity-level data, we can explore questions regarding 
how microbiomes compare across individuals, 
species or environments, yet there is a paucity of 
studies comparing microbiomes across animal spe- 
cies. Ley et al. (2008) compared gut microbiomes of 
60 different mammal species and discovered that 
both host phylogeny and diet appear to influence 
the bacterial communities. Roughly half of the 
samples they observed were from animals in the 
wild, but from populations that did not overlap in 
their geographic range, and the other half were from 
captive zoo animals. Yildirim et al. (2010) compared 
the gut microbiome of three different primate 
species from wild troops whose home ranges over- 
lapped. While the authors found highly significant 
differences in the gut bacteria among the primate 
species, it is difficult to completely disentangle the 
role of social interactions (for example, grooming, 
mating, aggressive interactions) and close physical 
contact. It is possible that the similarity of the gut 
communities of primate individuals is, in part, due 
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to closeness and physical contact within a troop 
rather than phylogenetic host differences. Thus, 
many wild animal species may not directly overlap 
in their habitat use or level of social contact with 
other species, and these factors may confound our 
understanding of species-specific differences of 
their microbiomes. 

As the skin microbiome is in more constant 
contact with the external environment than the 
gut microbiome, we might expect this to obscure 
patterns of phylogenetic host specificity. For 
example, we may expect the skin microbiome to 
consist of many exogenous bacteria originating from 
the environment (Belden and Harris, 2007), rather 
than bacteria that are host specific. Although the 
human skin microbiome and its spatial variability 
has been reasonably well studied (Costello et al., 
2009; Grice et al., 2009), patterns of microbial 
assemblages living on the skin of other animals is 
largely unknown. Without fur or feathers, amphi- 
bians provide an excellent model system for 
studying skin-associated microbial communities. 
Pre-metamorphic amphibians in pond habitats 
provide a rare example of multispecies communities 
that co-exist temporally and spatially in aquatic 
habitats. We can be confident that all larval indivi- 
duals present in a given pond during a breeding 
season originated there and are exposed to the same 
microbial inocula over the course of their larval 
development. If there are no differences in skin 
bacterial communities across species that co-habit 
in a pond, but there are differences across pond 
sites, then we can argue that the composition of 
skin-associated bacterial communities is most 
strongly influenced by the pond environment. 
Conversely, if we observe that the skin communities 
are unique to different amphibian species and that 
those patterns are consistent across multiple sites, 
then we can assert that the skin bacterial commu- 
nities are demonstrating evidence for host specifi- 
city. Our results indicate strong support for the latter 
case, skin bacteria of amphibians appear to be 
strongly species specific across different sites, even 
when species co-exist in the same pond (Table 2). 

The number of unique skin bacterial phyla varied 
between 10 and 18 phyla depending upon amphi- 
bian species. Using similar techniques, Costello 
et al. (2009) found the human skin microbiome is 
comprised of 18 different bacterial phyla. Therefore, 
the phylum-level bacterial richness found on the 
skin of amphibians appears high and on par with 
that of human skin. Previous studies have reported a 
diversity of bacteria from amphibian skin using 
different methods. Using a PCR-based fingerprinting 
technique (denaturing gradient gel electrophoresis) of 
16S rRNA fragments, Lauer et al. (2007) reported up 
to 19 unique bacteria from eastern red-backed 
salamanders [Plethodon cinereus), representing four 
different phyla. Likewise, Woodhams et al. (2007b) 
used a culture-based approach and found 40 unique 
bacterial isolates, representing three different phyla, 
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from 70 frog [Rana muscosa) individuals (mean of 
eight isolates per frog individual). As in other 
environments, culture-based approaches likely cap- 
ture only a small fraction of the bacterial diversity 
present on amphibian skin (Pace, 1997). While 
culture-based approaches allow one to determine 
the specific attributes of individual bacterial taxa, 
such as their ability to limit fungal growth (Lauer 
et al., 2007, 2008), combining molecular with 
culture-based approaches is clearly important for 
understanding the role of bacterial communities on 
their amphibian hosts. 

All of the phyla reported from previous studies of 
bacteria on amphibian skin that employed different 
methods, including Firmicutes, Actinobacteria, 
Proteobacteria and Bacteroidetes (Culp et al., 2007; 
Lauer et al., 2007, 2008; Woodhams et al., 2007a; 
Lam et ah, 2010), were also found in the present 
study and comprised some of the more abundant 
groups observed. Given that the barcoded pyrose- 
quencing method allowed us to observe 14 addi- 
tional phyla, we can examine community-level 
patterns and composition to a far greater degree. 
Closer examination of the composition of the 
bacterial communities revealed that, for several 
amphibians sampled in this study, their skin 
bacterial communities were dominated by a few 
phylotypes (Figure 3). For example, on tiger sala- 
manders and chorus frogs from four different 
pond sites, between 42-97% of the bacterial com- 
munity was dominated by members of the genus 
Curvibacter, a Betaproteobacteria in the family 
Comamonadaceae. This was by far the most abun- 
dant phylotype across all species in all sites. This 
genus and its close relatives are commonly found in 
a wide variety of natural aquatic systems (Shaw 
et al., 2008; Mueller-Spitz et al., 2009; Hutalle- 
Schmelzer et al., 2010) and they appear to be 
primarily freshwater organisms with a wide range 
of temperature tolerance (9-40 °C) and pH tolerance 
(5.5-8.5) for growth (Ding and Yokota, 2004). Due 
to the lack of comparable studies, we do not know 
if members of this genus are commonly associated 
with other species of amphibians found at other 
locations. Another dominant phylotype was identi- 
fied as being related to Ideonella [Betaproteo- 
bacteria] and was particularly abundant on western 
chorus frogs. The sequence matches that of a 
symbiont found on the epithelium of Hydra, a small 
freshwater cnidarian (Fraune et ah, 2010). A few of 
the other more abundant phylotypes found on the 
amphibians also match sequences of bacterial 
strains that are previously known as symbionts 
living with hosts. These include a strain of Sphin- 
gomonas [Alphaproteobacteria] reported to be a 
symbiont with plant hopper insects (Tang et al., 
2010), a strain of Hydrogenophilus [Betaproteo- 
bacteria] reported from the human oral microbiome 
(Dewhirst et al., 2010) and a strain of Caulobacter 
[Alphaproteobacteria] reported from the human gut 
(Frank et al., 2007). 



Overall, several bacterial types that are related to 
those found in freshwater systems (for example, 
Nevskia; Pladdies et al., 2004) are abundant on 
amphibian skin, as well as several that are related to 
known symbionts of other organisms. While we 
used sterile water rinsing of amphibians specifically 
to remove transient microbes before swab sampling 
(sensu Culp et al., 2007), it is possible that not all of 
the bacteria obtained from the swab are actual 
epibiotic symbionts of the amphibians. However, 
since we found that amphibian species inhabiting 
the same pond harbored distinct bacterial commu- 
nities, it is likely that most of the taxa were epibiotic 
and not simply transient microbes from the pond 
environment. Furthermore, we do not know whether 
transient microbes potentially interact with the 
epibiotic community on amphibians or simply exist 
temporarily in proximity to the skin microbiome. 
The mechanism by which amphibian species 
acquire unique skin communities is not known. 
Sources of microbial inocula may include vertical 
transmission (from parent to offspring), horizontal 
transmission (between individuals) and a variety of 
environmental sources (water, sediment, algae, in- 
vertebrates, organic matter, etc.). While some terres- 
trial species of amphibians exhibit parental care of 
egg nests and young (Banning et al., 2008), the 
aquatic species studied in this system do not display 
parental care; thus, vertical transmission seems 
unlikely in this system. However, Fraune et al. 
(2010) provide evidence that freshwater cnidarian 
[Hydra] embryos obtain maternally derived peptides 
that drive the assembly of the epithelial bacterial 
community, indicating that it may be possible for 
organisms without extended parental care to receive 
vertical transmission of peptides that select for the 
type of community that resides on the offspring. In 
any case, the assembly of microbial communities 
that live on amphibians is an intriguing area that can 
be examined with both field observations and 
experimental approaches. 

While Curvibacter is clearly the most abundant 
bacteria observed on amphibians in this data set, 
the presence and relative abundance of other 
phylotypes largely account for the species-specific 
patterns observed (Figure 3). Pond sites were 
not significant in determining the differences among 
individual amphibians (Table 2); however, we can 
see that there is some variation across sites; for 
example, chorus frogs from four different sites 
have variation in their communities (Figure 3). 
Three of the Colorado sites occur on the Front 
Range and one is on the Western Slope (~200 miles 
away), yet several dominant phylotypes are still 
found in both regions (for example, Curvibacter, a 
Sphingobacteriales phylotype, a Bacteroidaceae 
phylotype and Sphingomonas; see Figure 3). 
It is possible to study larger-scale regional varia- 
tion (for example, continental scale) for some 
wider ranging species such as leopard frogs, 
to see if dominant groups shift at different spatial 
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scales; however, these patterns are completely 
unknown. 

Skin-associated bacteria have the potential to 
influence susceptibility to infectious diseases that 
affect amphibians. Several studies have cultured 
bacteria from various amphibians and demonstrated 
that certain isolates can inhibit the growth of fungi, 
including Bd (Harris et al, 2006, 2009b; Lauer et al, 
2007, 2008; Woodhams et al, 2007b). Harris et al. 
(2009a) used an smti-Bd isolate [Janthinobacterium 
lividum) in a bioaugmentation experiment and 
demonstrated success in reducing mortality rates 
due to Bd for susceptible amphibians in the 
laboratory. The same bacterial strain and other 
anti-5d bacteria were found to be more abundant 
naturally on endangered mountain yellow-legged 
frog [Rana sierrae) populations that have co-existed 
with Bd relative to populations that have been 
extirpated due to Bd (Woodhams et al, 2007b; Lam 
et al, 2010). Just as some populations of amphibians 
differ with respect to their tolerance of Bd, different 
species of amphibians are known to vary in their 
ability to tolerate and resist infection with Bd 
(Garner et al, 2006; Kilpatrick et al, 2010). 
Our finding that amphibians have species-specific 
bacterial communities on their skin further suggests 
that these skin bacterial communities have the 
potential to influence species differences in suscept- 
ibility to Bd. Amphibians that harbor more anti- 
fungal bacterial species are likely to be better at 
tolerating or resisting infection with Bd. Broadly, 
how microbiotas influence the susceptibility to 
infectious disease is a next frontier with many 
potential applications relevant to human medicine 
(Stecher et al, 2010) and wildlife conservation 
(Harris et al, 2009b; Lam et al, 2010). 
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